Stability of inverse bicontinuous cubic phases in lipid-water mixtures 
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We investigate the stability of seven inverse bicontinuous cubic phases (G, D, P, C(P), S, I—WP, 
F — RD) in lipid-water mixtures based on a curvature model of membranes. Lipid monolayers are 
described by parallel surfaces to triply periodic minimal surfaces. The phase behavior is determined 
by the distribution of the Gaussian curvature on the minimal surface and the porosity of each 
structure. Only G, D and P are found to be stable, and to coexist along a triple line. The 
calculated phase diagram agrees very well with experimental results for 2:1 lauric acid/DLPC. 



Most of the many mesomorphic phases formed by lipid- 
water mixtures are of the inverse type, with the self- 
assembled lipid monolayers curving towards the aqueous 
regions H. In inverse bicontinuous cubic phases, a single 
lipid bilayer extends throughout the whole sample, divid- 
ing it into two percolating water labyrinths. Until now, 
the structures G, D and P have been identified. The only 
known lipid-water system in which G, D and P coexist 
is 2:1 lauric acid/DLPC and water M. The property of 
cubic lipid bilayer phases to divide space into interwoven 
polar and apolar compartments is utilized for biological 
function, e.g. in mitochondria and the endoplasmic retic- 
ulum ||. Recently, they have also been used as artificial 
matrices which enable membrane proteins like bacteri- 
orhodopsin to crystallize in a three-dimensional array M . 

It was shown by Luzzati and coworkers H that the 
mid-surfaces of the lipid bilayers are very close to cubic 
minimal surfaces, which have zero mean curvature every- 
where. These surfaces occur in lipid-water systems due 
to the local symmetry of the lipid bilayer, which implies 
that the surface should curve to both sides in the same 
way. However, it is well known U that many more cubic 
minimal surfaces exist than the structures G, D and P 
identified in lipid-water mixtures. What might be the 
reason that these other phases have not been observed so 
far? Helfrich and Rennschuh M have argued, based on 
the curvature model for fluid membranes [p|, that struc- 
tures with a narrow distribution of Gaussian curvature 
over the minimal mid-surface should be most favorable. 
However, the relevant data was known to them only for 
G, D and P, which are degenerate in this respect due 
to the existence of a Bonnet transformation. Recently, 
we obtained numerical representations for a large num- 
ber of cubic minimal surfaces in the framework of a sim- 
ple Ginzburg-Landau model |1. In this Letter we use 
this data to investigate seven inverse bicontinuous cubic 
phases (G, D, P, C(P), S, I - WP, F - RD). For an 
illustration of the G and S surfaces, see Fig. [jL We find 
that the width of the different distributions of Gaussian 
curvature is indeed smallest for G, D and P and larger 
for all other structures considered. This proves for the 



first time why only G, D and P should be observed ex- 
perimentally. 

Our detailed investigation of the stability of bicontinu- 
ous cubic phases shows that the existence of the Bonnet 
transformation implies that with increasing water con- 
centration, G, D and P coexist along a triple line. We 
show that this sequence is determined by an universal ge- 
ometrical quantity, the topology index, and that higher 
order terms in the curvature energy and thermal mem- 
brane undulations do not lift this degeneracy. Further- 
more, we calculate phase diagrams as function of concen- 
tration and temperature, which are in good agreement 
with an intermediate temperature range of the experi- 
mental phase diagram for 2:1 lauric acid/DLPC and wa- 
ter§. 

The main contributions to the free energy of an inverse 
bicontinuous cubic phase are the curvature energy of the 
lipid monolayers and the stretching energy of the hydro- 
carbon chains pXJ ] . The curvature energy is described by 
the Canham- Helfrich Hamiltonian |§ 



U= I dA {2k (H - c Q ) 2 + R K) 



(1) 



where H and K are mean and Gaussian curvature, re- 
spectively, and the integral runs over the neutral surface 
of the monolayer. The model parameters are sponta- 
neous curvature cq, bending rigidity K and saddle-splay 
modulus k. In the following, positive mean curvature cor- 
responds to curvature towards the aqueous regions. The 
stretching energy can be taken to be harmonic about the 
average chain extension I. In order to attain optimal free 
energy, the neutral surfaces of the two monolayers should 
realize at the same time constant mean curvature cq and 
constant distance I to the minimal mid-surface. Since this 
is not possible geometrically, the system will always be 
frustrated. The relative distribution of the frustration on 
bending and stretching energies depends on I. In a simple 
polymeric model (which certainly overestimates the chain 
compressibility of lipid monolayers), bending rigidity and 
stretching modulus scale as I 3 and I -1 , respectively [ pd[ , 
and the ratio of the two energies scales as (l/a) 2 , where 
a is the lattice constant. For low water content, where 



a ~ 21, stretching is important and the bicontinuous cu- 
bic phases are suppressed by the hexagonal phase. In this 
Letter we consider the opposite limit, where the stretch- 
ing energy dominates and the layer thickness is nearly 
constant. We model the two monolayers as parallel sur- 
faces at distances ±1 away from the minimal mid-surface. 
Their geometrical properties are thus given by 



dA l =dA (l + Kl 2 ) ,H l 



-Kl 
1 + Kl 2 



,K l = 



K 



l + Kl 2 



(2) 



With these relations, the curvature energy of the lipid bi- 
layer can be expressed in terms of the geometrical prop- 
erties of the minimal mid-surface. 

Our analysis shows that each (inverse) bicontinuous 
cubic phase is characterized by two dimensionless num- 
bers, A and T, which are sufficient to describe its sta- 
bility and phase behavior. A is the width of the distri- 
bution of Gaussian curvature over the surface and will 
be discussed in more detail below. The topology index 
F can be motivated as follows. For a two-dimensional 
surface embedded in three-dimensional space, there ex- 
ist three motion-invariant and additive functionals of the 
shape: surface area A, Euler characteristic \ an d inte- 
grated mean curvature. For triply periodic minimal sur- 
faces, the integrated mean curvature vanishes. The val- 
ues for A and x depend on lattice constant and choice 
of unit cell, but one ratio can be constructed which does 
not. We define this universal quantitiy to be the topol- 

ogy index T = (A 3 /2ir\x\V 2 ) , where V is the volume 
of the unit cell. It describes the porosity of a structure 
(the larger its value, the smaller is the connectivity) and 
has been investigated extensively by Hyde and cowork- 
ers [g2[ . The gyroid G is the least porous of the known 
triply periodic minimal surfaces since it divides space into 
labyrinths with 3-vertices only. Then comes S which has 
both 3- and 4-vertices, and D, which has 4-vertices only. 
All other structures involve higher-coordinated vertices. 
In Table I we collect values for T for all cubic minimal 
surfaces considered in this Letter. The values for G, D, P 
and / — WP are exact since exact (Weierstrass) represen- 
tations are known. G, D and P are related by a Bonnet 
transformation, which changes the geometrical proper- 
ties globally (i.e. shape in embedding space, topology, 
space group and lattice constant) but not locally (i.e. 
metric and curvatures). The values of T for C(P), S and 
F — RD given in Table I follow from our numerical rep- 
resentations as isosurfaces of scalar fields $(r) which we 
obtained recently from a simple Ginzburg-Landau model 
for ternary amphiphilic systems M ; due to numerical lim- 
itations, they are not too accurate but can be expected 
to appear at the correct position in the overall sequence. 
For a given cubic minimal surface with surface area 
A* and Euler characteristic x P er conventional unit cell, 
the hydrocarbon volume fraction v of the corresponding 



inverse bicontinuous cubic phase is given by 

^/>7-'=-GH*(r 



(3) 



It turns out that for all seven structures considered the 
lattice constant in units of chain length is very well ap- 
proximated by a /I — 2 A* jv over the whole range of v 
(for v > 0.8, the bicontinuous structures begin to self- 
intersect). The average mean curvature of the monolay- 
ers as a function of v is 



(H^l = 



f dA l H l l _ (v/T) 2 

jdA 1 ~A-(v/ry 



(4) 



In the following, spontaneous curvature Cq is measured 
in units of l/l. Eq. (0) implies that a given spontaneous 
curvature is attained at the hydrocarbon volume fraction 
v = (4co/(l + Co)) 1 / 2 T. Hence if all structures consid- 
ered were stable, for increasing water concentration we 
would expect them to appear in a sequence of decreasing 
topology index, i.e. in the sequence G-S-D-I — WP 
- P etc. 

The stability of the different phases is determined by 
the free energy per unit volume. In the absence of ther- 
mal fluctuations, the free energy density can be written 
in dimensionless units as 
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where the integration extends over the minimal mid- 
surface in a conventional unit cell. Here we defined 
r = -k]2k and E(K*) = K*A*/8nx- The curvature 
model (Q) is well-defined for < r < 1. S(if*) varies 
over the surface, but can be shown to be invariant under a 
Bonnet transformation. Higher-order terms in the curva- 
ture energy of the monolayers would result in expressions 
involving higher powers of the combination "E,(K*)(v /T) 2 . 
The free energy density of the lamellar phase is fi, a = v. 

Eq. (^) can also be expressed in the form of an effec- 
tive curvature energy for the lipid bilayer (like Eq. M) 
for the monolayers). From an expansion for small chain 
length I, one obtains an infinite series in powers of Gaus- 
sian curvature K. In particular, the effective value for 
the saddle-splay modulus of the lipid bilayer is found to 
be Rbi — 2(k + 4cok). Since Ru increases linearly with the 
spontaneous curvature of the monolayers, the cubic bi- 
continuous phases are favored for cq > r/2, even though 
K of the monolayers is assumed to be negative. 

The first (bending) term in Eq. dq) measures the devi- 
ation of the monolayers' mean curvature from the spon- 
taneous curvature. It follows from Eq. (J3) that a narrow 
distribution of Gaussian curvature K over the minimal 
mid-surface translates into a narrow distribution of mean 



curvature H l over the monolayer, which is favorable in 
terms of the bending term, as pointed out by Helfrich 
and Rennschuh M. For a more detailed analysis, we note 
that the bending term is an area average of a complicated 
function of its Gaussian curvature K, which is distributed 
non-uniformly over the surface. We proceed by first mea- 
suring the distribution f(K) of Gaussian curvature over 
the minimal surface. Even in the case of a known Weier- 
strass representation, the function f(K) is analytically 
not tractable and has to be determined numerically in the 
form of a histogram {if,, /j}. Its knowledge allows to re- 
place the area integral in Eq. (||) by a (one-dimensional) 
integral over Gaussian curvature. We obtained {Ki, /;} 
for G, D, P and / — WP from their Weierstrass represen- 
tations and for C(P) , S and F — RD from our numerical 
representations. A short description of the procedure in 
both cases and plots of the histograms obtained are given 
in Ref. M . In order to quantify the width of the differ- 
ent distributions, we calculate A = ((K — (K)) 2 )/{K) 2 
where the average is over the surface and normalized. 
Like the topology index T, A is independent of scaling 
and choice of unit cell. Moreover, it can be shown that 
it is also invariant under Bonnet transformations. The 
values for A for the seven structures considered are given 
in Tab. I. We find that it is the smallest for G, D and 
P and higher for all other structures. Thus we conclude 
that although for given cq the different structures have 
different optimal values for v according to their topology 
index, only G, D and P will have small frustration due 
to the relatively small variation of mean curvature over 
their monolayers. 

Since the quantity S(i*T*) in Eq. (pj) is invariant under 
a Bonnet transformation, the free energy densities for G, 
D and P obey a scaling form vg((v/T) 2 ), with a uni- 
versal function g{x) which can be read off from Eq. (|5|). 
For small r, this function has g(0) = 1, a minimum at 
x ~ 4co/(l + Co), and diverges at x ~ 4. An analyti- 
cal Maxwell construction for this scaling function shows 
that the Bonnet-related structures form a triple line, such 
that the coexisting values of v scale with T. Moreover the 
triple line includes an excess water phase at v — 0. The 
coexistence with an excess phase is known as emulsifica- 
tion failure from surfactant systems pi : if the swelling 
with water reaches the point where the structure adapts 
best to the given spontaneous curvature, any additional 
water is just expelled into an excess phase. 

In addition to the bending energy, thermal fluctua- 
tions have to be taken into account for the calculation 
of phase diagrams. In the case of the lamellar phase, the 
membrane undulations give rise to a steric repulsion p3[ , 
which in our case reads 



Is 



32c 2 



k B T 



(6) 



(1-v) 2 

where c^ = 0.106 [OJ. In the bicontinuous cubic phases 
the fluctuations of the mid-plane lead to a renormaliza- 



tion of k and R |15|]. However, since k multiplies the 
average mean curvature squared, which vanishes for the 
mid-plane, only the renormalization of R at length scale 
£ has to be taken into account. We identify the typical 
length scale of a cubic structure with {K)^ 1 / 2 . This im- 
plies £/l = 2T/v, so that r in Eq. (||) gets renormalized, 
with 



TR 
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(7) 



Therefore the scaling form of the free energy remains 
unchanged both by higher order terms in the curvature 
energy and by thermal undulations, although the scaling 
function itself changes. 

In Fig. 0, we plot the free energy densities for cq — 1/6, 
r = and kJIzbT = 10. As predicted from our general re- 
sults above, we find that G, D and P are indeed the only 
stable cubic phases and that they become stable along 
a triple line, which is the lower straight line. The same 
line also marks the Maxwell construction for the emulsi- 
fication failure, the coexistence of P with an excess water 
phase. These curves for different values of Cq allow the 
construction of phase diagrams for given r and n/hsT as 
shown in Fig. |3j. Note that without any further interac- 
tions, D is only stable along a line in our calculation. We 
therefore do not show two-phase coexistence regions be- 
tween the three bicontinuous cubic phases in Fig. H, but 
only the intersections of the free-energy curves. With in- 
creasing Co, the regions of stability shift to higher values 
of v as predicted by Eq. (|j) . 

The existence of the triple line implies that subtle dif- 
ferences in the free energies, which arise from a more de- 
tailed treatment of stretching contributions, of thermal 
fluctuations, or of van der Waals or electrostatic interac- 
tions, can be expected to destroy the degeneracy between 
G, D and P. This explains qualitatively why all three 
cubic phases are rarely observed in the same system. 

We can now compare our results with the experimen- 
tal phase diagram for 2:1 lauric acid/DLPC and water, 
the only known experimental system in which G, D and 
P coexist. Since the head group region is much smaller 
than the tail region of the lipid, we can identify the hy- 
drocarbon volume fraction v with the lipid concentra- 
tion. Further, a linear relation between Co and T can be 
assumed, since it is well established in non-ionic surfac- 
tant systems. In both phase diagrams, the sequence L 
-G-D-P-W + P is observed with increasing water 
concentration pw ~ 1 — v. For r = 0.1, the tempera- 
ture range from T = 35°C to T — 50°G in the phase 
diagram of 2:1 lauric acid/DLPC and water corresponds 
roughly to the range 0.18 to 0.26 in spontaneous curva- 
ture Co- Experimentally, the hexagonal phase dominates 
at higher temperatures and large hydrocarbon volume 
fraction, and fluidity is lost at lower temperatures (with 
the main transition at 30°G). An extrapolation of the 



linear fit shows that the balanced temperature is located 
near 0°C, i.e. well below the temperature of the main 
transition. 
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FIG. 1. One of the two monolayers in the G (left) and S 
(right) structures, for hydrocarbon volume fraction v — 0.5. 
Both inverse bicontinuous cubic phases have spacegroup laid. 
However, only G is found to be stable both in experiments and 
theory. 
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FIG. 2. Free energy densities as a function of hydrocarbon 
volume fraction v for cq — 1/6, r = and KJksT = 10. The 
phases can be identified from the sequence of curves at the 
top, which corresponds to F - RD, C(P), I - WP, S, P, D, 
G and L from left to right. Only G, D and P are stable; they 
form a triple line (see inset) which includes a coexistence with 
excess water (lower dashed line). The second dashed line is a 
Maxwell construction between G and the lamellar phase. 
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FIG. 3. Numerically calculated phase diagram for r = 0.1 
and K/fegT = 10 as a function of water volume fraction 
pw — 1 — V and spontaneous curvature Co. The hatched 
regions indicate two-phase coexistences. 
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D 


P 


I-WP 


S 


F-RD 


C(P) 


A 


0.219 


0.219 


0.219 


0.483 


0.587 


0.650 


0.842 


r 


0.767 


0.750 


0.716 


0.743 


0.795 


0.654 


0.656 



TABLE I. Universal quantities, which characterize the ge- 
ometry of minimal surfaces. A is the width of the distribution 
of Gaussian curvature and determines the overall stability of 
bicontinuous cubic phases. The values for A are the same for 
G, D and P due to a Bonnet transformation. The topology 
index V describes the porosity and determines the position in 
the phase diagram. 



